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Abstract 

In this paper, we analyze various critical transmitting/sensing ranges for connectivity and coverage 
in three-dimensional sensor networks. As in other large-scale complex systems, many global parameters 
of sensor networks undergo phase transitions: For a given property of the network, there is a critical 
threshold, corresponding to the minimum amount of the communication effort or power expenditure by 
individual nodes, above (resp. below) which the property exists with high (resp. a low) probability. 
For sensor networks, properties of interest include simple and multiple degrees of connectivity/coverage. 
First, we investigate the network topology according to the region of deployment, the number of deployed 
sensors and their transmitting/sensing ranges. More specifically, we consider the following problems: 
Assume that n nodes, each capable of sensing events within a radius of r, are randomly and uniformly 
distributed in a 3-dimcnsional region 1Z of volume V, how large must the sensing range Rsense be to 
ensure a given degree of coverage of the region to monitor? For a given transmission range Rtransi what 
is the minimum (resp. maximum) degree of the network? What is then the typical hop-diameter of 
the underlying network? Next, we show how these results affect algorithmic aspects of the network by 
designing specific distributed protocols for sensor networks. 

Keywords 

Sensor networks, ad hoc networks; coverage, connectivity; hop-diameter; minimum/maximum degrees; 
transmitting/sensing ranges; analytical methods; energy consumption; topology control. 

I. Introduction 

Advances in micro-electro-mechanical systems technology, wireless communications and 
digital electronics have enabled the development of multi-functional sensor nodes. Sensor 
nodes are small miniaturized devices which consist of sensing, data processing and com- 
municating components [1], [18], [48]. These inexpensive tiny sensors can be embedded 
or scattered onto target environments in order to monitor useful informations in many 
situations. We categorize their applications into law enforcement, environment, health, 
home and other commercial areas. Moreover, it is possible to expand this classification 
with more categories including space exploration [28] and undersea monitoring [51]. We 
refer here to the survey paper [1, Section 2] for an extensive list of possible applications 
of sensor networks. These critical applications introduce the fundamental requirement of 
sensing coverage that does not exist in traditional ad hoc networks. In a sensing-covered 
network, every point of a targeted geographic region must be within the sensing range of 
at least one sensor. 



February 1, 2008 



DRAFT 



3 

In general, networked sensors are very large systems, comprised of a vast number of 
homogeneous miniaturized devices that cooperate to achieve a sensing task. Each node 
is equipped with a set of sensing monitors (light, pressure, humidity, temperature, etc.) 
accordingly to their designated tasks and uses radio transmitters in order to communicate. 
We refer, for example, to the web-sites [7], [11], [44], [58] of some research institutions 
dedicated to the study and development of these networks. Typically, a sensor node is 
able to sense events within a given radius (the sensing range). Similarly, any pair of sensors 
are able to communicate if they are within a distance less than their transmitting range of 
each other. Wireless sensor networks are usually multihop networks as opposed to wireless 
LAN environments. 

A commonly encountered model of sensor network is defined by a pair n and TZ where n 
homogeneous sensor nodes are randomly thrown in a given region 1Z of volume V = \1Z\, 
uniformly and independently. This typical modeling assumption is commonly used by 
many researchers [8], [12], [13], [22], [24], [25], [34], [49], [54], [55]. In particular, the initial 
placement of the nodes is assumed to be random when the sensors are distributed over 
a region from a moving vehicle such as an airplane. As opposed to traditional ad hoc 
networks, a sensor network is normally composed of nodes whose number can be several 
orders of magnitude higher than the nodes in an ad hoc network. These sensor nodes are 
often deployed inside a phenomenon. Therefore, the positions of the nodes need not be 
engineered or pre-determined. Note that many existing results focused on planar networks 
[12], [38], [49] while three-dimensional settings reflect more accurately real-life situations 
[25] , [28] . Three-dimensional networks arise for instance in building networks where nodes 
are located in different floors. 

The aim of this paper is two-fold: (a) to study the role that randomness plays in sensor 
networks and (b) to investigate the design and analysis of appropriate protocols for these 
networks. These issues are motivated by the following simple reasons. On first hand, the 
random placement of the nodes allows rapid deployment in inaccessible terrains. On the 
other hand, this implies that the network must have self-organizing capabilities. Therefore, 
we start deriving analytical expressions to characterize the topological properties of sensor 
networks. Next, we discuss how to use these fundamental characteristics in order to design 
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and analyze fundamental algorithms such as those arising frequently in classical distributed 
systems. To name a few, these algorithms include the broadcasting and gossiping protocols 
[4], [5], the leader election algorithm [36] and distributed code assignments protocols [27]. 

II. Related work 

Study of random plane networks goes back to Gilbert [22] in the early 60's. In compar- 
ison to the well-known and well-studied random Bernoulli graphs 1 G(n,p) [9], [17], [30], 
only few papers considered the probabilistic modeling of the communication graph prop- 
erties of wireless sensor or ad hoc networks and the theory of random geometric graphs 
(RGG) is still in development with many problems left open (cf. [47]). In contrary to 
the Bernoulli graphs, in random geometric graphs the probability of edge occurrences are 
not independent making them more difficult to study. In RGG, a set of n points are 
scattered, following a given distribution, in the deployment region 1Z. Penrose studied the 
longest edge of the Euclidean minimum spanning trees (MST) and the longest nearest- 
neighbor link [45]. The same author established also that if the region of deployment is 
a (i-dimensional cube then the graph of communication becomes /c-connected as soon as 
its minimum degree reaches k [46]. In a slightly different setting, Gupta and Kumar [24] 
stated that when n nodes are distributed uniformly in the disk of unit area and their trans- 
mitting range is set to r = ^ lnw +^( TO 2 then the resulting network is connected with high 
probability 2 if and only if c(n) — > oo. Their results are obtained making use of the theory 
of continuum percolation [39], which is also used in [16] to investigate the connectivity of 
hybrid ad hoc networks. 

As far as we know, the critical transmitting range and the critical coverage range have 
been investigated first in [49] for the case when the nodes are distributed in a square 
according to a Poisson point process of fixed intensity. For the line of a given length, it 
has been studied in [50]. In [55], Shakkottai et al. considered connectivity, coverage and 
hop-diameter of a particular sensor grid. 

In this paper, we consider a model similar to those of [49], [50], [52], but for three- 

x E.g. there are more than 800 references in [9]. 

2 An event £ n is said to occur with high probability or asymptotically almost surely (a.a.s. for short) if its 
probability tends to 1 as n — > oo. Formal definitions will be shortly given. 
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dimensional settings. Furthermore, our results are illustrated with two fundamental dis- 
tributed protocols intended for the neighborhood discovery and the code assignment prob- 
lem [6], [27]. In particular, our results permit to do the average-case based analysis of the 
code assignment problem which has been investigated empirically in [6, Section 4]. 

The remainder of this paper is organized as follows. Section 3 presents formal de- 
scriptions of the models used throughout the paper. Section 4 offers results about the 
relationships between the sensing (resp. transmitting) range, the number of nodes and 
the volume of the region to be monitored. In particular, we show how to quantify the 
minimum and maximum degrees of the network. We also show how to compute the hop- 
diameter of the underlying graph (defined as the maximum number of hops between any 
pair of sensor nodes), when the transmitting range is slightly greater than the one required 
to have a connected network. In Section 5, we present general algorithmic schemes for 
distributed protocols related to our settings. In particular, we consider how to design a 
polylogarithmic protocol to allow the nodes to discover their neighborhoods asymptoti- 
cally almost surely. Next, we turn on the design and analysis of a protocol to assign codes 
in such random networks. The code assignment problem consists to color the nodes of a 
graph in such a way that any two adjacent nodes are assigned two different colors. 

III. Preliminaries and models 

To analyze the topology of sensor networks, three fundamental models are needed: (a) a 
model for the spatial node distribution, (b) a model for the wireless channel of communica- 
tion between the nodes and (c) a model to represent the region monitored by a particular 
sensor. Throughout this paper, nodes are randomly deployed in a subset 1Z of IR 3 following 
a uniform distribution. To model the wireless transmission between the nodes, a radio link 
model is assumed in which each node has a certain communication range, denoted R TB ans- 
Two nodes are able to communicate if they are within the transmitting range denoted 
R-trans of each other. Only bidirectional links are considered. Clearly, one imposes a graph 
structure by declaring any two of the stations that lie within a given transmission range 
to be connected by an edge. The resulting graph is called the reachability graph and is 
denoted G(p, R TB ans) where p — y represents the (expected) density of the nodes per unit 
volume. The network coverage is defined as follows. For a node located at a point p, 
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its monitored region is represented by a sphere, the sensing sphere, of radius R Sekse and 
centered at p. A region 7Z is said covered if every point in 7Z is at distance at most R SENSE 
from at least a node. We assume that every node has the same sensing range R SENSE and 
the same transmitting range R Tb .ans- Given the increasing interest in sensor networks, our 
purpose is to design a solution where the nodes can maintain both sensing coverage and 
network connectivity. 

Throughout this paper, V always represents the volume of the region 7Z to monitor and 
we allow both n, V — > oo, in such a way that n/V — > p (dense networks). Concretely, 
p = y represents the expected number of nodes per unit volume. 

Notations 

• Degrees of the reachability graph. The degree of a node v represents the number 
of its neighbors in the graph G(p, R TrA ns) an d is denoted d v . We consider here questions 
such as: What is the required value of the transmitting range R Tr ans to have a reachability 
graph with a given minimum (resp. maximum) degree 5 (resp. A)? 

• Diameter. The hop distance between two nodes u and v is defined as the length of 
the shortest path (with respect to the number of hops) between them. The diameter (or 
hop- diameter) of a graph is the maximum of minimum hop distances between any two 
pair of nodes. Under the same hypothesis as above, what is the typical diameter of the 
reachability graph of the random sensor networks? 

• Degrees of coverage. Define any convex region of 1R 3 as having a degree of coverage 
k (i.e., being k- covered) if every point of the considered region is covered by at least k 
nodes. Given the volume V of 7Z, what is the required value of the sensing range R SENSE to 
achieve a specified coverage degree k, k > 0? 

Throughout this paper, we will use standard mathematical notations [15] concerning 
the asymptotic behavior of functions. We have 

• f(n) = 0(g(n)) if there exists a constant c and a value n such that f{n) < cg{n) for 
any n > n . 

• f(n) = Q(g(n)) if there exists constants c\ and C2 and a value no such that c\g{n) < 
f(n) < C2g{n) for any n > uq. 
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• f(n) = o(g{n)) or f(n) < g(n) if ^ -> as n -> oo. 

• /(n) ~ if and only if lim n ^oo = 1. Equivalently, f(n) ~ <7(n) if and only if 

f(n) = g(n) + o(g(n)). 

• We recall that an event £ n (depending on the value of n) is said to occur asymptotically 
almost surely (a.a.s.) if its probability tends to 1 as n — > oo. We also say £ n occurs with 
high probability (w.h.p.) as n — > oo. 

IV. Fundamental characteristics of random sensor networks 

To generate the network, the n sensor nodes are distributed in the fixed region of volume 
V . If p — y tends to some constant (i.e., p = 0(1)), this process is well approximated by 
a Poisson point process of finite intensity p (see [26]). Note that this assumption is well 
suited for both theoretical point of view and in practice since by a suitable rescaling [39] 
all the properties obtained in this paper can be reformulated. For instance, any realization 
G(p, Rtrans) using a transmission range R TBANS and with expected number of nodes p per 
unit volume coincides with another realization G(p', Rtrans') with a transmission range 
r' Trans and intensity p' provided that p' = (Rtrans/R TR ans) 3 P- This model is also well suited 
for faulty nodes since if each node is independently faulty with some probability p, the 
properties obtained here remain valid with n replaced by p x n. 

A. Connectivity regime and minimum transmission range 

As a warm-up, we begin with a natural question that often arises concerning the mini- 
mum value of Rtrans required to achieve connectivity [24]. We have the following theorem : 

Theorem 1: Let 1Z be a bounded and connected set of M 3 of volume V . Suppose that 
n nodes are placed in 1Z according to the uniform distribution and assume that p — y — 
0(1). The network formed by adding edge between nodes of distance at most 

R-Thans = 

S/3(lnn + cu(n))/47rp (1) 

is connected if and only if lim n ^ 00 cj(n) = oo. 

Proof: By assumption, the distribution of the nodes can be approximated (see [26, 
Section 1.7]) by a Poisson point process of intensity p — y which has the following property. 
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The probability pe(i) that a randomly chosen node % has t neighbors is given by 

Ptii) = (4/3P 7 LJ ' ex P (" 4 /3p-R? RANS ) • (2) 



Thus, if Rtrans = \/3(lnn + uj{n)) /Anp, as n — > oo no sensor node is isolated with proba- 
bility close to 

J](l-po(i)) = (l-e- 4/3 ^—j =h-i_- J =ex P UlnM j j 

= expf-e-^+of^^^y (3) 



n 

Therefore, there is a.a.s. no isolated nodes if and only if uj(n) tends to infinity with n. To 
prove that the graph is also connected, we argue as in [24]. More precisely, we make use 
of results from continuum percolation [39]. In percolation theory, nodes are distributed 
with Poisson intensity p and as in our model, two nodes are connected iff the distance 
between them is less than r{p). Denote by G Poisson (p ; R(p)) the resulting infinite graph. 
Theorem 6.3 of Meester and Roy [39] states that almost surely G Poisson (p, R(p)) has at 
most one infinite-order component for p > 0. Thus, almost surely the origin (the node 
distribution is conditioned on the origin having a node) in G Poisson (p ; R(p)) lies in either an 
infinite-order component or is isolated. Our problem can be approximated by regarding 
that process as the restriction to 1Z of the Poisson process of intensity p on IR 3 . Denote 
by G^ isson (p, R(p)) this restriction of G Poisson (p ? R(p)) to 1Z. By the above observation 
(see also [24, Section 3] or [16]), the probability that G^ oisson (p, R(p)) is disconnected is 
asymptotically the same as the probability that it has at least one isolated node. For large 
n, the difference between the G^ isson (p, R(p)) and our model can be neglected. Thus, to 
ensure connectivity it suffices that there is no isolated vertex. By (3), this is only achieved 



upon setting R TrA ns = \/3(lnn + u{n)) /Anp, with u(n) — > oo. ■ 
Remark. It is important to note here that Theorem 1 concerns all connected regions 1Z 
of R 3 of bounded volume (V = \1Z\ < oo). 

B. Coverage and minimum sensing range 

To study coverage properties, we need some results from integral geometry [53], [42]. 
The following lemma is due to Miles [43]. 
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Lemma 2: Let 7Z be a bounded set of M 3 of volume V . Let P be a point process of M 3 of 
intensity p. Suppose that each point x of the point process V monitors a sphere centered at 
x and of radius r. Define N(p) as the number of these spheres containing p e ffi 3 . For any 
set 5 C M 3 , denote respectively by N~(S) = inf p&s iV(p) and by N + (S) = sup pgiS iV(p). 
Then, the following holds 



lim sup 



and 



lim sup 

V->oo r>0 

with 



and 



P 



[N-(K) > j] -exp ( - 3pvr 2 /32 j(j + 1) ^ <p(j + h 4vrpr 3 /3)) 



0, (jGN) 
(4) 



P 



[JV+(ft) <€]- exp (-f(i,p,V) (1-^-1, 4vrpr 3 /3)) ) | ^ , G N*) (5) 



p, V) = p (4 + 3/8(tt 2 + 1Q)(£ - 1) + 3tt 2 /32(£ - !)(£ - 2)) V 



ip(x,y) = e 




(6) 



(7) 



Lemma 2 tells us about the limiting probability distribution of the number of spheres 
covering a point p of V . For a given sensing range R SeN se, to ensure the total coverage of 
the region 1Z (that is every point of TZ is within the sensing range of at least one sensor), 
it suffices to have N~(1Z) > 1. The following result establishes the relation between R SENSE , 
V and n : 

Theorem 3: Let 1Z be a bounded set of 1R 3 of volume V. Assume that n sensors, each of 
sensing range R SENSE , are distributed uniformly and independently at random in 71. Suppose 



that p 



y is constant. Let 



R 



Sense 



s/ (Inn + lnlnn + uj(n)) V 



Ann 



(8) 



With probability tending to 1 as n — > oo (and V = 0(n)), every point of 71 is monitored 
by at least one sensor if and only if u{n) — > oo. 

Proof: We have to find the value of R SENSE (in terms of n and V) such that as n — > oo (or 
V — > oo), the probability P [TV - (7?.) > 0] — > 1. By formula (4) in lemma 2, this probability 
is well approximated by exp (— 3p7r 2 /16 V <p(2, 47rr 3 /3)) (we used j = 1 in the formula). 
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lnn+ln lnn+iv(n) 



Therefore, it suffices to have nip (2, 47rr 3 /3) <C 1. By setting R SENSE = ^3 , . ( 
we obtain 

nip(2, 4npr 3 /3) = nip(2, Inn + lnlnn + uj{n)) 

= n exp (— Inn — In Inn — u){n)) 

( 1 2 

x I f + Inn + In Inn + cu(n) + — (Inn + In Inn + u)(n)) 



(9) 



exp (u(n)) 

Hence, P [N-(K) > 0] -> 1 if and only if w(n) -> oo. ■ 

Remark. Formulae (1) and (8) indicate that the radius R SEKSE required to achieve a sensing- 
covered network is greater than the transmission range R TbA ns required to have a connected 
network. For randomly deployed sensor nodes, these results show that if the transmis- 
sion/sensing radius can be "compared", connectedness arises slightly before sensing cov- 
erage. It is also important to note that the arguments used to prove these Theorems show 
that both connectivity and coverage properties are subject to abrupt change known as 
phase transition phenomena. From an engineering standpoint, these formulae are crucial: 
For instance, (1) indicates that there is a critical transmission effort required to each in- 
dividual node to ensure with high probability that any pair of nodes in the network can 
communicate with each other through multihop paths. We refer here to the paper of Goel 
et al. [23] and to Krishnamachari et al. [32] for results about phase transition phenom- 
ena in wireless networks and random geometric graphs. In particular, the authors of [23] 
proved that monotone properties for random geometric graphs have sharp thresholds (we 
refer to Friedgut and Kalai [20] for the classification of critical thresholds phenomena in 
random graphs). 

C. Degrees and coverage 

In various networking problems, multiple-paths between any pair of two nodes are of 
importance. On one hand, the existence of multiple independent paths plays crucial 
role (e.g. in flooding or in gossiping protocols). On the other hand, in dense networks 
where each node has a great number of neighbors, the number of interferences makes the 
scheduling of communications difficult. 
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Similarly, different applications may require different degrees of sensing coverage. While 
some protocols require that every location in the considered region be monitored by one 
node, other applications need significantly higher degrees of coverage. In these directions, 
one is first interested in 1-coverage (which is solved by Theorem 3) but also in /c-coverage 
for several values of k, in particular for unbounded value of k, viz. k = k{n) ^> 1. Our aim 
in this paragraph is to address issues related to the probabilistic relationships between the 
transmission/sensing ranges, the number of nodes and the size of the region of deployment. 
More precisely, we are interested in values of the sensing range R SENSE such that the event 
"each point p of the considered region 1Z is covered by at least k spheres of radius R SENSE " 
occurs asymptotically almost surely. Similarly, we are interested in several degrees of 
connectedness depending on the transmission range. 

Theorem ^: Let 1Z be a bounded subset of M 3 of volume V . Assume that n sensors 
are deployed uniformly and independently at random in 71 and p = y = 0(1). Then as 
n — > oo, the following holds : 

(i) Given any fixed integer £ > 0, if the sensing range R SENSE satisfies 4/3 7rpR 3 ENSE = 
Inn + £ In Inn + uj{n) then 1Z is a.a.s. ^-covered if and only if 1 <§; oj{n) <C In Inn. 

(ii) Given any function c(n) satisfying 1 C c(n) C j^il. jf 4/3 n p r 3 ^ = In n+c(n) In In n 
then each point of 1Z is a.a.s. monitored by at least ~ c(n) spheres and at most ~ e Inn 
spheres. 

(Hi) For any constant real number £, £ > 0, if 4/3 7rpRg ENSE = (1+^) Inn then the number 
N(p) of sensors covering and monitoring each point pofTZ satisfies a.a.s. 



where W (z) denotes the branch 3 of the Lambert W function W(z) satisfying —1 < W(z) 
and W-i(z) is the branch of the same function satisfying W(z) < —1. 
(iv) Given any function c(n) satisfying 1 <^ c(n) <^ if 4/3 7rpR 3 ENSE = c(n) Inn then 
each point p of 1Z IS db.db. s. covered by N(p) ~ c(n) Inn. 

Proof: The proof of (i) is very similar to the one of Theorem 3 and is therefore 
omitted. The proofs of (ii), (iii) and (iv) rely on extremely precise analysis of the 

3 See the Appendix for details on the Lambert W function [14]. 




£ Inn 



£ Inn 



+ o (Inn) , 



(10) 
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truncated gamma function p present in equations (4) and (5) above. More specifically, the 
gamma function is denned as 



oo 

tj.X-1. 



T(x) = / e~H x - x dt. 







The incomplete gamma functions arise from the integral above by decomposing it into two 
integrals, the first one from to y and the second from y to oo : 



7 (x,y) = I" e'H^dt, 9*e(x)>0, 
Jo 

/•oo 



T(x,y) = / e~ l t x - l dt, |argy|<7r. (11) 
Jy 

Specific cases are obtained when x is an integer. For x G N, we have 

T(x + l,y) = x\e-« ( £ £ ) , (12) 
\fe=o ' / 

Thus, for any natural number y the function ip described in (7) can be expressed as 

T(x + l,y) 

^ V)= T(x+l) ■ (13) 
Uniform asymptotic expansions for the so-called incomplete gamma function are now 
required to cope with the value of tp present in (4), resp. (5). These expansions were 
derived by Temme [56] (see also [21] for a survey on incomplete gamma functions). By 
using the integral representation of the truncated gamma, viz. 

T(r v) e xc t>(P) r c + io ° i / \ dt 

1 {X ' y) = - / e xm 0<c<8 (14) 

r(x) 2m J c _ lO0 (3-t> u<c <^ 

where 

(j)(t) =t-l-\nt and (3 = - , (15) 
the author of [56] arrives at an asymptotic representation of the form 

m 2 y\i2j 

e~2 xr > / 1 l\ e~2 xr i I Ci(rj) 



In (16), the erfc is the complementary error function defined by 

2 f x 

erfc(x) = 1 - — / e^ 2 rft, (17) 
V 71 " io 
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the coefficients Ciijf) can be computed by a power series expansion (we refer to [56] for 
details) and 

In particular, as x — > oo if rj is bounded the remainder terms of (16) is exponentially small. 
Therefore, we have 



T(y) 2 V\2J ^/2^\P-l vj \V2nx 

We quote here that computations of the functions q, involved in (16), as well as the 
methods to get them are detailed in [56], [57]. Now using (16) and letting the value of j 
(resp. £) vary, we can make the values 

exp (-3 P 7r 2 /32i(j + l)V<p(j + 1, 47rpr 3 /3)) (20) 

in (4) and its counterpart 

exp(-p (4 + 3/8(tt 2 + 1Q)(£ - 1) + 3tt 2 /32(£ - 1)(£ - 2)) V (l - <p(£ - 1, 47rpr 3 /3))) 

(21) 

in (5) very close to 1. In fact, suppose for instance that 4/37rpR| ENSE = Inn + c(n) In Inn 
as given in the hypothesis of (ii). Therefore, ip(j + 1, 4/3 up R| ENSE ) = (p(J + l,mn + 
c(n) In Inn). Since the expression inside the exponential in (20) must tend to 0, the value 
of <f(j + 1, Inn + c(n) In Inn) must behave like 

ip(j + 1, In n + c(n) In Inn) = -— [ , (a > 0) . (22) 

v v ' ' fn lnn a 



erfc(x) = (l + O (1 ] ), *-oo (23) 



Using the fact that, 



with the help of formula (19), we then have to solve (asymptotically) with respect to j 

1 e -\ii+i)v 2 ( i \ / n\\ 

1 + , , (24) 



fn\nn a y/2n(j + 1) \P - 1 / V \3 

which implies that (5 > 1. Equivalently, we have to find j satisfying 



J 3„l n „« V 1 + (jJJ=-H-2 y + 1) "J' <25) 
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Using (18) and taking the logarithm of the previous equation, we obtain 

(j + 1) (P - 1 - ln/3) - - ln j = Inn + a lnlnn - ln(/3 - 1) + O (1) . (26) 

Consequently, we then have : 

1 V ^/lnlnn\ ^ /ln/3\ ^/lnA 

Mnn + O — — )+0 i-f )+0 i-f . (27) 



p-i-\npj v. P J \ P J \ P 

Clearly, the equation (27) above can be solved asymptotically (with respect to j). To this 
end, we need results concerning the asymptotic behavior of the Lambert W function [15]. 
Whenever z approaches 0~, W-i has a complete asymptotic series expansion provided by 
the Lagrange inversion theorem (see [14] and [15]) which starts as 

W^z) = In (-*) - In (- In (-*)) + O i^^f^) ' (* "> ' z < °) " ( 28 ) 
Substituting /3 by ln n+c( "! ! ln ln n and using standard analysis, it yields for j in (27) : 



j'+i 

(n) lnlnn 



(l + o(l))=w(n) + o(w(n)) , (29) 



yj/ / cj(n) lnlnn \ 

— 1 I e(lriTi+iij(n) lnlnn) J 

where we used the function W-i(z) and the asymptotic formula (28) for W-i(z) as z — > 0~. 

Now, we turn on the upper-bound - given by (ii) - of the number of spheres monitoring 
each point of 1Z. We argue as above. Fix a constant a > 0, not necessarily the same as 
in (22). As for formula (22), the expression inside the exponential in (21) can be made 
sufficiently close to if we find £ such that 

fn (1 - <p(i - 1,4/3 Tip rLJ) = , (30) 

which implies that 

,n<9 1 4/^n R 3 ^ r(l,4/3 7T P Rf EKSE ) 1 

^-l,4/37r P R SENSE ) = ^ = 1-^-^-. (31) 

Instead of (23), we now have the following asymptotic expansion : 

erfc(-x) = 2 - + O ( 1— | , x -> oo . (32) 
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Then, by setting 4/3 7r p Rg ENSE = Inn + c(n) In Inn, we find 

i = e Inn + o (Inn) (33) 

as stated for the upper-bound given by the statement of property (ii) in the Theorem. 

The prove (iii), we first remark that for any function u(n) 3> 1, the functions W and 

W_i involved above verify 

win) u(n) , . 
^~ — — r-~w(n), n^oo. (34) 

Denoting by j\ the at least number of spheres monitoring each point of 1Z when setting 
the sensing range to 4/3 7rpR| ENSE = (l + £) Inn, instead of (27) we have to solve 

3i= ( l + ln( l ^ J J lnn + 0(lnlnn) , (35) 

Equation (35) can be solved asymptotically (w.r.t. ji) and this time we find 

ji — 7 v + (l nn ) • (36) 



e (1+*) / 

Similarly, denote by ji the at most number of spheres covering each point of 1Z, we can 
argue as (33) to find j 2 : 

J 2 = t ^ + o(lnn). (37) 



e {l+i) / 

Combining (36) and (37), we just have the property (iii) as stated. Similarly, one can 
prove the property (iv) of the Theorem. ■ 

Observe that the previous results reflect also the degrees of the nodes of the reachability 
graph simply by replacing the sensing range R SeN se with the transmitting range R TrA ns and 
we have the following : 

Corollary 5: Let 1Z be a subset of M 3 such that \1Z\ = V < oo and assume that n sensor 
nodes are deployed in 1Z following a uniform distribution. Suppose that y is a constant 
p. Denote by 5 (resp. A) the minimum (resp. maximum) degree of a given network and 
for any node v, d v denotes the degree of v. The network formed by adding edge between 
nodes of distance at most R Tr ans have the following properties asymptotically almost surely 
as n — > oo : 
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(i) For any fixed integer £ > 0, if the transmitting range R TBANS satisfies 4/3 7rpR^, BANS = 
Inn + £ In Inn + u(n) then the minimum (resp. maximum) degree 5 (resp. A) of the 
reachability graph satisfies 5 > £ + 1 (resp. A < elnn). 

(ii) For any function c(n) s. t. 1 «C c(n) <C if 4/3 7rp R 3 RANS = In n + c(n) In Inn then 
each node v of the network has a degree d v comprised between c(n) and e In n. 

fmj For any constant real number £ (£ > 0), if 4/3 7rp R^ RANS = (1 + Q Inn then the degree 
d v of any node v of the network satisfies 

£ Inn . £ Inn , . , 

+ o(\nn)<d v < -. v + o(\nn), (38) 



('iv^ Given any function c(n) s. t. 1 c(n) if 4/37rpR^ RANS = c(n)lnn then each 

node v of the network has a degree d v ~ c(n) In n. 

Proof: For each one of the statements (i), (ii), (iii) and (iv), by substituting R SEKSE 
in Theorem 4 with R TrA ns, each geographical point p of 1Z lies inside the corresponding 
number of "spheres of communications" . In particular, the points representing the centers 
of the sensors are asymptotically inside the same number of spheres. Therefore, the proof 
of the statements (i), (ii), (iii) and (iv) in the Corollary follows the previous proof of 
Theorem 4. ■ 

D. Hop-diameter 

In this paragraph, we consider a setting slightly different from the above. For sake of 
simplicity, we consider a cubic region 1Z of the Euclidean space R 3 . In what follows, the 
transmission range of the nodes is set to a certain value R TrA ns such that 

A = lim 4/3 ™ Rtrans , A e (0, oo] . (39) 
n^oo V Inn 

For A defined by (39) and borrowing terms from Bernoulli random graphs G(n,p) [9], [17], 
[30], with respect to the connectivity property, the regime is referred to as : 

• the subcritical regime if A < 1, 

• the critical regime if A = 1 and 

• the supercritical regime if A > 1 . 

For transmission ranges of the form R TRANS = ^/S/An (1 + £) lnn/n, that is the connectivity 
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regime is supercritical, we have the following result for the order of the diameter of the 
network: 

Theorem 6: Suppose that n sensor nodes are randomly deployed in a cubic region of 
volume V of 1R 3 according to the uniform distribution. If their common transmission 
range is set to R TrA ns = \J ^^Ir ^" with t > 11/5, then the diameter D of the network 
satisfies : 

r ~ " (40) 



lim P 

n— >oo 



Proof: 



D < 12 i 3 



6(l + f) Inn 



1. 




Split the cube of volume V into k 3 sub- 
cubes, IZi, 7ZkA of equal volume. 
Each of these sub-cubes has side \fV / k. 
Choose k such that a sphere of radius 
Rtrans can entirely fit inside a sub-cube 
IZi (cf. figures). 



That is (under the hypothesis of the Theorem), k = ^ 



R (T7tt-\ — • F° r sake of 

6 Inn 

simplicity, let us suppose that this value of k is an integer. By (38) of Corollary 4, with 
high probability there is 6 (Inn) nodes inside the sphere of radius V 1 ^ 3 /2k. Consider now 
two adjacent sub-cubes as depicted in the figure below : 

The lens-shaped region, here denoted 
Li, represents the intersection of two 
spheres of radius R TKANS and whose cen- 
ters are at distance R TrA n.s- A bit of 
trigonometry shows that the volume of 




such intersection is given by 



5ttRt 



12 



According to the uniform distribution, there is no node inside each lens of volume \L\ 



5ttR t 



12 



with probability 



j^il 
V 



5(l + £) lnn\ n 



16n 



(41) 
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Since each sub-cube has at most 6 lenses similar to L 1; none of these lenses is empty with 
probability at least 

/ / 5 (1 + £) lnn\ n \ 6k3 ( ( ( 5 (1 + 1) lnn\ \ \ 6fc3 

/ / 5 (1 +£) \nn\ \ m"'- 

5 —)) 

£ exp (-(TT^ x jUPm) ■ (42) 

Therefore, if £ > 11/5 each lens similar to L x is non-empty a.a.s. Hence, from two adjacent 
sub-cubes TZi and IZj, communications between any node uGTJj and any node v e 72.^- 
need at most 6 hops as shown by the figure below 




By simple counting arguments, the 
proof of the Theorem is now complete. 



Remark. We conjecture that in the supercritical regime for connectivity, that is for 
transmission ranges of the form R TRANS = {/3/47T (1 + £) \nn/n with < £ < 11/5 the 
diameter of the network is (in probability) of order of magnitude 0(n 1 / 3 / lnn 1//3 ). In 
contrary, it would be much more difficult to capture the diameter in the just critical 
regime, viz. if R TRAm = (Inn + uj{n))/n with u{n) ^> 1. 

V. Distributed protocols 

In this section, we consider some distributed protocols which are built on the top of 
the previous results concerning the main characteristics of random sensor networks. As 
emphasized by the fundamental papers [3], [4], [5], to cite only a few, the number of 
nodes n, the diameter D and the maximum degree A of the networks can play crucial role 
when designing distributed protocols for radio networks. In fact, the executions of these 
algorithms are often measured in terms of n, D and A. 
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A. A model for protocols 

A commonly encountered model for distributed protocols is described briefly in this 
paragraph. We refer to [3], [5] for detailed descriptions of this model. A distributed 
protocol for multihop networks is a protocol executed at each node in the network in the 
following way : 

• The time of execution is considered to be slotted and are subdivided into time slots or 
rounds. 

• In each round, a node acts either as transmitter or as receiver. A node u receives 
a message m sent by one of its neighbors in a given round if and only if (a) it acts as a 
receiver and (b) exactly one of its neighbors acts as a transmitter. If two or more neighbors 
of u are sending at the same time-slot the node u does not receive nothing. That is, nodes 
are unable to distinguish between collision and the lack of message. 

• The nodes are assumed to be distinguishable, that is each node has an unique identifier, 
ID for short, ranging from 1 to n. In our settings, we assume also that the nodes are aware 
of their number n as well as the volume V of the region of interest. 

Remark. In what follows and without loss of generality, we assume that the nodes can 
transmit messages up to a distance of order 0{yj^^-). We note that if the (common) 
transmission range of the nodes is such that R TrA ns 3 ^ hl „ V , the results presented in the 
following paragraphs can be easily extended using the same global ideas. 

B. A simple protocol for neighborhood discovery 

The first distributed protocol which will be discussed is a protocol called ExchangeID. 
This algorithm allows each node u to discover the set of its neighbors, denoted T(u). This 
protocol appears to be useful since as already stressed the nodes are deployed in a random 
fashion and therefore, they do not have any a priori knowledge of their respective neigh- 
bors. The neighborhood discovery is done using a simple randomized greedy algorithm 
whose pseudo-code is given in the following : 
(0) Protocol ExchangeID (n, V, R TbA ns) 



(1) Begin 

(2) Com 
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(3) Then set A := -I x lnn/W (-^) and G t := 2 exp (2^0 (-^)) 5 

(4) For i := 1 to [Qlnn 2 ] Do 

(5) With probability each node u sends a message containing its own ID ; 

(6) EndFor 

(7) End. 



Theorem 7: Suppose that n sensor nodes are randomly deployed in a region of vol- 
ume V following a uniform distribution. If their transmission range satisfies R TbA ns > 
[ln(n)+u(n)]V, with u(n) > 1 but R TRANS = O (^/^P^), then after one exe " 
cution of ExchangeID (n, V, R Tr ans), with probability tending to 1 as n — > 00, every node 
has received correctly all the identities of all its neighbors. 

Proof: The proof of Theorem 7 relies on two facts, viz., (1) the main characteristics 
of the random Euclidean network and (2) the number of iterations T = C^ln(n) in the 
main loop of ExchangeID is sufficient for the nodes to send its ID at least once to all its 
neighbors. For the first point (1), we have seen that for any node v of the network, the 
degree of v (d v = |r(t>)|) satisfies w.h.p. (cf. Corollary 5 equation (38)) : 

dv< ! lnn x +o(\nn) . (43) 



e(l+e) 

Therefore, at the regime considered in the hypothesis of Theorem 7, the maximum degree 
of the graph is (with high probability) bounded by Qlnn (where q satisfies e.g. q = 
2 W (— i/e (1 + £)). Using this latter remark, let us complete the proof of our Theorem. 
For any distinct pair of adjacent nodes and any time slot t G [l, Ct In (n) 2 ] , define 
the random variable (r.v., for short) xfl- as follows: 

1 if and only if the node j does not receive the ID of % 
1,1 - 1 at time t G [l, C t ln(n) 2 ] , (44) 

• otherwise . 



yW _ 



In other terms, the set 

{x^,t,j^t,te[l,C e \n(n) 2 ]} (45) 
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denotes a set of random variables that counts the number of "arcs" % — > j such that j has 
never received the ID of %. Denote by X the r.v. 

X = Y,*~i (46) 

where X^j = 1 iff xf\ - = 1 for all t G [l, Cg In (n) 2 ] . Now, we have the probability that 
i does not succeed to send its ID to j at time t: 

d 



P 



In (n) j ^ Inn I V Inn 



Therefore, considering the whole range [l, Cg In (n) 2 ] , after a bit of algebra we obtain 

P[X^. = 1] < (l 



-ce \ ln-(n) 2 C'i 



In (n) 

< exp (- In (n)e- ct C e ) (48) 

which bounds the probability that i has never sent its ID to j for all t G [l, ln(n) 2 C£]. 
By linearity of expectations and since by (38) the number of edges is of order 0(n\nn), 
we then have 

E[X] < O (nln (n) exp (- In (n) e- Ce C t )) ■ (49) 

Thus, E[X] <C 1 as n — > oo for a certain constant Q such that, say 

C< > 2e Ce . (50) 

Note that this constant can be computed for any given £ using, e.g. eg — 2 W {—i/e (1 + £)) 
and using the first moment method [2], one completes the proof of Theorem 7. ■ 

C. Coloring and codes assignment 

According to the rules of distributed protocols given above, when a node is transmitting, 
all the nodes within its transmitting range must be silent. Moreover in multihop environ- 
ments, collisions can also occur when two non-adjacent nodes are trying to transmit to a 
common neighbor [hidden collisions). To circumvent these problems, researchers use to 
assign orthogonal codes to the nodes, a problem equivalent to that of coloring the graphs 
associated to the physical network [6], [19], [27]. 
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To assign codes to the nodes of the network, let us consider the following simple and 
intuitive randomized protocol called AssignCode. Each vertex u has an initial list of 
colors (also referred as palette) of size d u + 1 = \T(u)\ + 1 and starts uncolored. We can 
assume that each node knows its neighbors (or at least a great part of them) by using 
the previous algorithm, viz. ExchangeID. Then, the protocol AssignCode proceeds in 
rounds. In each round, each uncolored vertex u, simultaneously and independently picks 
a color, say, c from its list. Next, the station u attempts to send this information to his 
neighborhood denoted by T(u). Trivially, this attempt succeeds iff there is no collision. 
Before attributing the color c definitely to u, its neighbors has to sent one by one a 
message of reception. Note that this can be done deterministically in time O(lnn) since u 
can attribute to its active neighbors in T(u) a predefined ranking ranging from 1 to 
Therefore, u sends a message of confirmation and its neighbors undergo an update of their 
proper palettes and of their active neighbors. Hence, at the end of such a round the new 
colored vertex u can quit the protocol. The details of ASSIGNCODE follows : 

( 0) Protocol AssiGNCoDE(n, V, R TsANS ) 
( 1) Begin 

( 2) Each vertex u has an initial palette of colors, say p{u) = {c±, c 2 , • • • , 



3) 
4) 
5) 
6) 
7) 




Then set A := - I x \nn/W 
For i := 1 to D t In (nf Do 




and choose e. g. D e : 



For each vertex u do 



8) 
9) 
10) 

11) 
12) 

13) 
14) 
15) 



• Pick a color c from p{u) ; 

• Send a message containing c with probability 



i 



If no collision Then 



Every station v in T(u) gets the message properly ; 

One by one (in order) every member of T(u) sends a message ; 

(* This step can be synchronized by always allowing A time slots. *) 



Endlf 




u sends a message of confirmation and goes to sleep ; 
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(16) 



every station in T(u) removes the color c from its palette ; 



(17) Endlf 

(18) EndFor 

(19) End. 

Theorem 8: Suppose that n sensor nodes are randomly deployed in a region of vol- 
ume V following a uniform distribution. If their transmission Tctngc satisfies R-trans ^ 
[\n(n)+u(n)]V, with u(n) > 1 but R TRANS = O (j/^£r) , then after one ex- 
ecution of AssiGNCODE(n, V, R T rans), with probability tending to 1 as n — > oo, every 
pair of nodes at distance at most R TRANS from each other have received two distinct codes 
(colors). 

Proof: Although more complicated, the proof of Theorem 8 is very similar to the one 
of Theorem 7. For any distinct node u, recall that T(u) represents the set of its neighbors 
and denote by p u the size of its current palette. Now, define the random variable Y u as 
follows: 



Denote by T„ the set of active neighbors of u at any given time t during the execution 
of the algorithm. Suppose that we are in such time slot t. Independently of its previous 
attempts, u remains uncolored with probability 



Y u ={ 



• 1 if and only if the node u, 

remains uncolored after the D e In n 2 



steps of AssignCode 



(51) 



• otherwise . 





There is at least a collision due to one neighbor v € V 



u 



(52) 
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Since Vt, |r^| < A and Vt>, 1 < \p v \ < A + 1, we have 



x ir (t) i 

' \ I / / 1 " 



^ i 1 - (at^j + r-^-A, 1 1 (53) 



i \ i / / 



* i 1 - <a7^) J + (AT^j x i 1 -I 1 - a) ' (54 > 
£ '-^aTSj - 1_ ^a - (55) 



Therefore, with probability at most 



Z^lnn 2 



1 \^ m " _ / D e lnn 2 \ 

'-eK) ^("^H (56) 

u remains uncolored during the whole algorithm. Thus, the expected number of uncolored 
vertices at the end of the protocol AssignCode is less than 

E[Y] = E[T U ] < n exp (-^^) ■ (57) 

Since by (38) we have 

A = A(Q<-2 (58) 

It is now easy to choose a constant .D^ such that 

12£ 

D »- Wo (- e/ e (! + ()) • (59) 

in order to have E[Y] <C 1 as n — > oo. After using the well known Markov's inequality (cf. 
[2]), the proof of our Theorem is now done. ■ 

VI. Conclusion 

The main purpose of this paper has been that of investigating the fundamental char- 
acteristics of a randomly deployed set of sensor nodes. First of all, with respect to the 
communications, the characteristics of interest include the diameter and the degrees (min- 
imum and maximum) of the reachability graph generated by the nodes. Next, taking 
the sensing range of the nodes as a parameter, several degrees of coverage have been 
characterized with rigorous proofs. 
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On the top of these typical behaviors of random sensor networks, two distributed mech- 
anisms are proposed. The first one concerns the dissemination of the identifiers of the 
nodes to their neighbors while the second protocol solves the code assignment problem. 
With high probability, both protocols are shown to achieve their tasks in polylogarithmic 
time slots. 

While performance evaluations of algorithms intended for networks have been employed 
for mainly simple and regular graphs, our results show that current mathematical tools 
are available to design and analyze protocols intended to sensor networks. 

Finally, we briefly point out a remark related to some real-life situations. A possible 
way of future investigations could be to search for similar results as those presented here 
whenever the communication radii of the nodes do not degrade immediately but rather in 
a continuous fashion. 
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Appendix 



The Lambert W function. In this paragraph, we give some properties of the function 
satisfying W(x)e w ^ = x. We remark here that the function W, in particular the princi- 
pal branch Wq, already plays a central key role when studying the random graph model 
G(n, m), i.e., the random graph built with n vertices and m edges which is the "enumer- 
ative counterpart" of the G(n,p) random graph model (see, e.g., the "giant paper" [29]). 
In fact, — Wq(— x) is the exponential generating function that enumerates Cayley's rooted 
trees [10] and we have 



We plot in Figure 2 the two real branches of the Lambert W function considered in 
this paper. This function has been recognized as solutions of many problems in various 
fields of mathematics, physics and engineering as emphasized in [14]. The Lambert W 
is considered as a special function of mathematics on its own and its computation has 
been implemented in mathematical software as Maple. Figure 2 represents the two real 
branches of the Lambert W function. It is shown that the two branches meet at point 
M = (— 1/e, —1). As an example, if £ — \ in (10) of Theorem 4, each point of the area X is 
covered, with high probability, by at least j ~ .1520088850 Inn disks. We have the Figure 
1 depicting the functions £ — > — i/W-i e ^ +t ^ j and £ — > — £/W ^— e( ^ +t) j involved in 
Theorem 4. 




(60) 



i=i 
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Fig. 1 

Plots of -i/W-i (-^^jy) and 
-W(-^y). 



Fig. 2 

The branches Wo (dashed line) and 
W-i (solid line) of the Lambert W 

FUNCTION. 
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